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PREFACE 


This  technical  note  is  the  third  in  a  series  on  the  utilization  of 
the  computer  to  assist  mutagenicity  testing  at  Letterman  Army 
Institute  of  Research.  In  this  report  we  detail  three  functions: 
searching  of  a  specific  file  to  summarize  the  data;  calculating  the 
mean  mutation  frequency  with  its  variance  expressed  as  percentages; 
and  using  a  statistical  program  for  analysis  of  the  data.  Using  the 
program  in  the  Drosophila  melanogaster  Sex-Linked  Recessive  Lethal 
Assay  provides  an  easy  manner  to  view  the  data  and  perform  statistical 
analyses. 
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Computer-Assisted  Statistical  Analysis  III. 

The  Drosophi la  melanogaster  Sex-Linked 
Recessive  Lethal  Assay 

The  Drosophila  melanogaster  Sex-Linked  Recessive  Lethal  Assay 
requires  the  testing  of  8000  to  10,000  X-chromosomes  for  each  compound 
and  concurrent  negative  control  in  addition  to  500  X-chromosomes  for  a 
positive  control.  A  mass  of  data  must  be  processed.  Therefore,  two 
computer  programs  were  developed  to  summarize  the  results.  By 
totalling  the  number  of  lethal  mutations  and  non-lethal  mutations  for 
each  brood  of  each  replicate  yielding  the  total  number  of  lethal 
mutations  and  number  of  tests,  this  first  program  calculates  the 
percentage  mean  mutation  frequency  and  its  variance.  The  second 
program  Derforms  the  statistical  analysis  using  the  BMDP  statistical 
software™  package  by  testing  the  null  hypothesis  if  the  number  of 
lethal  mutations  resulting  from  the  compound  is  non-significantly  from 
different  from  the  number  resulting  from  the  negative  control. 

DESCRIPTION  OF  THE  PROGRAM 

The  program  DRUNSUM.FR  (Appendix  A)  searches  the  data  file  and 
summarizes  the  data  accumulated  in  the  program  DROSTOXDATA.FR  (1)  for 
a  specified  subgroup  of  a  given  run  number  selected  by  the  user.  This 
program  uses  the  subroutines  RDDROSDATA.FR  and  STRING. FR  (Appendix  B 
and  C)  to  read  the  data  sequentially  and  accumulates  the  values 
necessary  for  the  number  of  lethal,  non-lethal,  and  total  tests 
performed  for  each  brood,  and  grand  total  for  the  specified  run 
number.  The  program  also  calculates  the  mean  mutation  frequency  and 
its  variance  and  displays  them  as  percentages.  The  detection  of  three 
error  conditions  is  provided  to  see  (a)  if  more  than  25  records  are 
present  for  the  specified  subgroup  within  the  run;  (b)  if  the  file 
name  under  which  the  data  from  the  DR0ST0XDATA  are  stored  is  present; 
(c)  if  the  specified  subgroup  is  present. 

The  program  DRUNSUM.FR  (2)  requests  the  following  information  from 
the  user  (Figure  1); 

•  The  name  of  the  file  to  be  searched. 

•  The  run  number  to  be  searched. 

•The  group  designator  of  the  male  code  number  (this  specifies 
if  the  male  was  from  the  control,  positive,  or  test  compound). 

*PU 

The  second  program  is  part  of  the  BMDP  statistical  software  (3). 
Instructions  for  its  use  are  given  in  this  statistical  package.  This 
package  analyzes  the  data  in  a  2  x  2  table  and  two  statistics  are 
used,  the  Yates  corrected  chi-squares  and  Fisher’s  Exact  test.  The 
Yates  correction  is  intended  to  improve  the  approximation  to  the 
chi-square  distribution.  The  Fisher's  Exact  test  is  computed  when  the 
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minimum  value  is  less  than  20.  Both  1-tail  and  2-tail  probabilities 
are  tested.  However  the  2-tail  probability  for  Fisher's  Exact  test 
should  be  used. 

The  program  requests  the  following  information  non  the  user: 

•Enter  the  number  of  lethal,  lethal  mutations  and  non-lethal 
.nutations  resulting  from  the  negative  control. 

•Enter  the  number  of  lethal  mutations  and  the  number  of 
non-lethal  mutations  resulting  from  the  concurrent  test  compound. 

The  computer  will  then  type  the  values  and  probat  1 1 i ties  associated 
-ith  the  Yates  corrected  chi-square  and  the  Fisher's  Exact  test. 

COMMENT 

By  utilizing  these  programs  and  the  subroutines  presented  in  this 
report,  the  raw  data  may  be  summarized  and  statistical  analysis 
performed . 

CONCLUSION 


None 


LfcCOMMENDATION 

None 
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SAMPLE  RUN  OK  DR UN 


)  DR  UN 

The  computer  will  type: 

NAME  OF  FILE  TO  BE  SEARCHED:  DR OS DATA 


The  computer  will  type: 

••••DROSOPHILA  RECORD  SEARCH  AND  SUMMARY**** 
PLEASE  ENTER  RUN  NUMBER  TO  BE  SUMMARIZED:  24 

The  computer  will  type: 

ENTER  GROUP  TO  BE  SEARCHED  FOR: 

(T 1 ,  T2,  PI  or  Cl  ETC.  ALLOWED):  £1 

The  computer  will  type: 


SUMMARY  TOTALS 

RUN  NUMBER:  24 

GROUP  TYPE:  Cl 

BROOD 

DATA 

LETHALS:  0 

1 

0 

1 

NONLETHALS:  596 

524 

547 

547 

GRAND  TOTAL  LETHALS:  2 
GRAND  TOTAL  NONLETHALS: 

2214 

TOTAL  RECORDS  PROCESSED:  25 

TOTAL  TESTS  REPRESENTED:  2216 

MEAN  MUTATION  RATE:  .08000%,  WITH  VAR.  .07667% 

The  computer  will  then  display  the  user  may  then  return  to 

the  program  for  the  next  file  to  be  searched  or  log  off. 


Figure  1.  Sample  Run  of  DRUN 
(Underlined  entries  are  sample  user  inputs) 
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APPENDIX  A 


C  ***  PROGRAM  NAME:  DRUNSUK.FR 

C  ***  WRITTEN  BY:  WARREN  W.  JEDERESRG 

C  ***  PURPOSE:  TO  SUMMARIZE  THE  DATA  FROM  SPECIFIED 

C  ***  RUNS  IN  DROSDATA. 

C  ***  THE  SUBROUTINE  RDDROSDATA  IS  USED  TO  READ  EACH  RECORD 
C  ***  VARIABLE  DEFINITIONS 
C  N  “NUMBER  OF  RECORDS 

C  FBL  “FIRST  BROOD  LETHAL3 

C  FBNL=FIEST  BROOD  NONLETHALS 

C  STL  “SECOND  BROOD  LETKALS 

C  STNL=SECOND  BROOD  NONLETHALS 

C  TTL  “THIRD  BROOD  LETHALS 

C  TTNL=THIRD  BROOD  NONLETHALS 

C  FTL  “FOURTH  BROOD  LETHALS 

C  FTHL-FOURTH  BROOD  NONLETHALS 

C 

C  GTT  =  GRAND  TOTAL  TESTS 

C  GTL  “GRAND  TOTAL  LETHALS 

C  GTNL“GRAND  TOTAL  NONLETHALS 

C  TITO  “TOTAL  MUTATION  RATE 

C  VI KR“ ACCUMULATED  MUTATION  RATES  SQUARED 

C  NNR  “DESIGNATED  RUN  NUMBER  (USER  INPUT  AT  RUN  TIME) 

C  NS  “SUBGROUP  DESIGNATOR  (USER  INPUT  AT  RUN  TIME) 

C  RMR  “MEAN  MUTATION  RATE 

C 

C  VARIABLES  MENTIONED  IN  THE  COMMON  STATEMENT  ARE  DEFINED  IN 
C  THE  FORTRAN  V  PROGRAM  DROSTOXDATA. 

INTEGER  NR , MNUM , CNAM , BFF , BFL , 3FNL , 3SF , BSL , BSNL , BTF , BTL , 3TNL  f 
1 BLF ,  BLL ,  BLNL ,  TF ,  SL ,  ML , TL  f  TNL , TT ,  SFNAI-I 

DIMENSION  NS(2),SFNAM(20) 

REAL  MR 

COMMON/DREC/NR , MNUM( 6 ) , CNAM ( 6 ) , BFF , BFL , BFNL , BSF , BSL , BSNL , BTF , BTL , 
1 BTNL , BLF , BLL , BLNL , TF , SL , ML , TL , TNL , TT , MR 


v.ViV.y.AV.v. 
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APPENDIX  A 
(continued ) 


LOGICAL  EOF 

C  ***  SET  ACCUMULATORS  TO  0.00 
N=0 
FBL-0 
FBNL=0 
STL=0 
STNL=0 
TTL-0 
TTNL*0 
FTL=0 
FTNL=0 
GTT*0 
GTL*0 

gtnl=o 

TMR30 
VI  MR=*0 

C  ***  INITIALISE  EOF 
EOF3. FALSE. 

C  ***  GET  FILE  NAME  TO  CHECK  FOR  THE  DATA 
WRITE  (10,220) 

220  FORMAT  (  "'NAME  OF  FILE  TO  BE  SEARCHED:  ”,2) 

READ  (11,210)  SFHAM 
210  FORMAT  (2CA2) 

CALL  STRING  (SFNAM,20) 

OPEN  1 ,SFNAM,ATT3"SIB”,ERR3100 
READ  (1 ,5) 

5  FORMAT  (/) 

GO  TO  120 

100  WRITE  (10,6)  SFNAM 

6  FORMAT  (//,"****  ERROR  ****  ERROR  ****  ERROR  ****",/, 
14X,"FILE  ”,20A2,"  NOT  FOUND  .  .  .  .") 

STOP 

120  CONTINUE 

C  ***  PRINT  HEADER  AND  GET  RUN  NUMBER 
WRITE  (10,1) 

1  FORMAT  (/,1  OX, ''****  DROSOPHILA  RECORD  SEARCH  AND  SUMMARY 
1 10X ,  "PLEASE  ENTER  RUN  NUMBER  TO  BE  SUMMARISED:  ",i'J 


*#** 


ACCEPT  II NR 
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.APPENDIX  A 
(continued) 


C  *#*  GET  GROUP  DESIGNATOR 
99  CONTINUE 
WRITE  (10,3) 

3  FORMAT  (/.I OX, "ENTER  GROUP  DESIGNATOR  TO  EE  SEARCHED  FOR 
1 10X,“ (T1 ,T2,P1 ,0R  Cl  ETC.  ALLOWED) Z) 

READ  (11  ,4)(IJS(I),I-1  ,2) 

4  FORMAT  (2A1) 

IF  (HS(1).HE.MC  " . AND . NS ( 1 ) . NE . "P  ".AND. NS ( 1 ) .NE. "T  ")  GO  TO  99 
C  ***  READ  DATA  RECORDS 
140  CALL  RDDROSDATA  (EOF) 

IF  (EOF)  GO  TO  225 

IF  (NR.NE.NNR)  GO  TO  140 

IF  ( NS ( 1 ) .  NE . KNUM CO)  GO  TO  140 

IF  ( I1S  ( 2 ) .  NE .  HNUM(  2 ) )  GO  TO  HO 

N=N+1 

IF  (N.GT.25)  GO  TO  500 

FBL=FBL+BFL 

PBNL-FBNL+BFNL 

STL=STL+BSL 

STNL=STNL+BSNL 

TTL=TTL+BTL 

TTNL*TTNL*BTNL 

FTL=FTL+BLL 

FTNL=*FTNL+BLNL 

CTL'CTL+TL 

GTNL=GTNL*TNL 

GTT=GTT'fTT 

TMR=TMR+MR 

V1MR-V1MR*(KR**2) 

GO  TO  140 

C  ***  IF  N>25  ERROR 
500  CONTINUE 
WRITE  (10,7) 

7  FORMAT  (/,"***  ERROR  ***  ERROR  ***  ERROR  ***",/, 

1 4X ,  "MORE  THAN  25  RECORDS  FOUND 
STOP 
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APPENDIX  A 
(continued ) 


225  CONTINUE 

C  ***  IF  NO  SPECIFIED  RECORDS  FOUND 
IF( N.NE.O.O)  GO  TO  300 
WRITE  (10,9)  NNR.NS.SFNAM 
9  FORMAT  (/,"***  ERROR  ***  ERROR  *#*  ERROR  ***",/, 

1 4X , "RUN  NUMBER : ",2X, 14, 1 X, "GROUP:", IX, 2A1 ,/, 

24X , "NOT  FOUND  IN  " , 1 X ,20A2 , "  .  .  .") 

STOP 

300  CONTINUE 

C  ***  CALCULATE  MUTATION  RATE  AND  VARIANCE 
RHR-(TMR/H) 

VHR- ( VI  MR- ( ( TKR**2 . 0 ) /H ) ) / ( N- 1 ) 

C  ***  WRITS  OUT  TOTALS  • 

WRITE  (10,8)  HNR , NS , FBL , STL , TTL , FTL , FBRL , STNL, TTNL , FTNL , 

1 GTL, GTNL , N , GTT , RKR , VMR 

8  FORMAT  ( 20X , "  SUMMARY  TOTALS:”,/, 

1 1 OX , "RUT  NUMBER : " , 2X , 14 , / , 

210X, "GROUP  TYPE: " ,2X,2A1 ,/, 

320X ,  ’’BROOD  DATA:",/, 

410X, "LETHALS: " ,4(5X , 13) ,/ • 

57X , "NONLETHALS : " , 4 ( 5X , 13 ) , // , 

61 OX ."GRAND  TOTAL  LETHALS:" ,5X, 15,/, 

77X, "GRAND  TOTAL  NONLETHALS :",5X, 15,//, 

81  OX, "TOTAL  RECORDS  PR0CESSED:",5X,I2,//, 

910X, "TOTAL  TESTS  REPRESENTED:" ,5X, 15,/, 

1 1 0X , "KEAN  MUTATION  RATE:”,1X,F10.5,"S,  WITH  VAR.",1X,F10.5, 
CLOSE  1 
END 
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APPENDIX  B 


COMPILER  NOSTACK 
C  ***  PROGRAM  NAME: 

C  ***  WRITTEN  BY: 

C  ***  PURPOSE: 


RDDROSDATA.FR 

WARREN  JEDERBERG 

TO  READ  DATA  FOR  MAKING  REPORTS 


FOR  THE  SLRL-DROSOPHILA  ASSAY 


SUBROUTINE  RDDROSDATA  (EOF) 

COKMON/DREC/NR , MNUM( 6 ) , CHAM( 6 ) , BFF , BFL , BFNL , BSF , BSL , BSNL , 
1 BTF ,  BTL ,  BTNL ,  BLF ,  BLL ,  BLIIL ,TF,SL,ML,TL,  TNL ,  TT ,  MR 


INTEGER  NR , KHUK , CRAM , BFF , BFL , BFNL , BSF , BSL , BSNL , 3TFT2TL . 3TNL , 
1 BLF , BLL , BLNL , TF , SL , EL , TL , TN  L , TT 


REAL  MR 
LOGICAL  EOF 

READ  (1,100,  END*  120)  NR, KRUK.CKAM, BFF, BFL, BFNL, BSF, BSL, BSNL, 
1 BTF, BTL, BTNL, BLF, BLL, BLNL, TF,SL, ML, TL, TNL, TT, MR 

100  F0RKAT(I4,1X,6A1 ,1X,6A1 ,4(1X,3I2),U,I2,1X,3I2, 
12(1X,I3),1X,F6.2) 

RETURN 

120  EOF  -  .TRUE. 

RETURN 

END 


— -  TO  CONVERT  FILE  NAMES  TO  USE  IN  OPEN  STATEMENTS 
SUBROUTINE  STRING ( LINE, LLEN) 

COMPILER  STATIC 
DIMENSION  LINE  (LLEN) 

DO  100  1-1, LLEN 

IF  (LINE(l).EQ."  M)  LI1IE(I)=0 

IF  ( FLD( LINE( i) ,9,16) .EQ. FLD( "  " ,9, 1 6) )FLD(LINE(I) ,9, 1 6 )-0 

100  CONTINUE 
RETURN 


